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We report on the complete OPE series for the 6-gluon MHV and NMHV amplitudes in planar 
Af = 4 SYM theory. Namely, we provide a finite coupling prediction for all the terms in the expansion 
of these amplitudes around the collinear limit. These furnish a non-perturbative representation of 
the full amplitudes. 


In [T], building on [2], the so-called POPE (Pentagon 
Operator Product Expansion) program was put forward. 
Its goal is to put to full use the integrability of the chro¬ 
modynamic flux tube of planar Af = 4 SYM theory to 
study scattering amplitudes and polygonal Wilson loops 
at any value of the’t Hooft coupling g = VT/dTr. 

Since then, through a series of works [SHE], a simple 
physical and mathematical structure underpinning these 
objects has emerged. In this short letter, we collect all 
known information about hexagonal loops - i.e. 6-gluon 
scattering amplitudes - and present their complete non- 
perturbative expression as a POPE expansion. 

Within the POPE formalism, the renormalized hexag¬ 
onal Wilson loops mm are given by an infinite sum 



over all possible n-particle states of the colour flux tube, 
integrated over their n-dimensional phase space, and with 
Sn a symmetry factor for identical particles; it is equal 
to ni\n 2 \ ... Ufc! with n = ni + .. .rik and rij the num¬ 
ber of identical particles of a given kind. The POPE 
integrand in Q stands for the probability that the cor¬ 
responding state gets successively produced and annihi¬ 
lated on, respectively, the bottom and top pentagons in¬ 
side the hexagon. It can be conveniently decomposed as 
a product of three simpler building blocks as 

n = Ildyn X IIfF X Ilinat • (2) 


The product of the first two terms was spelled out ex¬ 
plicitly in [E|- The last factor, the so-called matrix part, 
deals with the 51/(4) R-symmetry degrees of freedom 
of the excitations, and was partially unveiled in for 
scalar excitations in the hexagon MHV Wilson loop. Be¬ 
low we quote this last missing factor for any state on any 
hexagon. The derivation of this matrix part factor and 
its generalization to any polygon will be given in m- 
The dynamical part neatly factorizes as 
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where r, a and cj) are the thee independent conformal 
cross ratios we denote in the OPE as time, space and an¬ 
gle, see [T]. The hexagon geometry enters only in these 
exponentials. Most importantly, we have the pentagon 
transition link P{ui\uj) between excitation i and excita¬ 
tion j and their corresponding flux tube measures giui). 

Next we have the form factor part given once more by 
a simple factorized product 


Hff = g 
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where r^ {rt) = 0,1,..., 4 is the R-charge carried by the 
bottom (top) pentagon [EHEj- For NMHV hexagons we 
have rt = 4 —r^, while for MHV and N^MHV amplitudes 
we have = r* = 0 and ri, = rt = 4, respectively. 
Therefore this factor is only needed to describe (any of 
the five independent) NMHV hexagons. 

Note that P, /i, h and the propagation charges implic¬ 
itly depend on the flavour of the excitation; e.g. if exci¬ 
tation j is a fermion and excitation j a scalar (j>, then 
in the notation of nni we have P{ui\uj) = Pijj 4 ,iui\uj), 
pi{ui) = jjL^iiiui), h{ui) = h^{ui), E(ui) = E,p{ui) etc. All 
form factors, transitions and measures are summarized 
in |12j while the dispersion relations can be found in |5]. 
For scalar and gluon excitations the contour of the in¬ 
tegration in Q is over the real axis. For fermions it is 
slightly more involved. As described in [3|, for these ex¬ 
citations it is often convenient to divide the integration 
contour into so-called large and small fermion contours. 

Finally we have the matrix part which encodes the 
SU{4) R-symmetry structure of the theory. It comes 
about from the contraction of R-charge indices between 
the bottom and top pentagons. The matrix part is a 
coupling independent rational function of (differences of) 
the rapidities {ui\. For generic number of excitations it 
does not exhibit any obvious factorization, in contrast 
with the two contributions described thus far. As an ex¬ 
ample, consider a three-scalar contribution in an NMHV 
hexagon with two units of R-charge at the bottom and 
similarly at the top. For such state the matrix part is 
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FIG. 1. The pattern of auxiliary rapidities entering the ma¬ 
trix part has a simple group theoretical interpretation. The 
three sets of rapidities can be identified with the three 

nodes of the SU (4) Dynkin diagram. The occupation numbers 
K\, 2,3 are fixed such that the overall state with scalars, 
Nij, fermions and anti-fermions belongs to the SU(4) rep¬ 
resentation { 1 , 4 , 6 , 4 , 1} for Tb = {0,1, 2, 3, 4} respectively. 
The cartoon depicted here is related to 0 by identifying the 
solid lines with /’s and the dashed lines with g’s. 


given by 

(7 + E s? - E (12 -b E s? - E 

i<j 

where {si, 52 , 53 } are the rapidities of the three scalars. 
The matrix part for states with more scalars or fermions 
typically lead to gigantic rational functions. Nonetheless, 
viewed from an appropriate angle, a neat factorization is 
actually there even for the matrix part! 

Indeed, it so happens that this matrix part itself ad¬ 
mits an integral representation over auxiliary rapidities 
with a very simple and totally factorized integrand. To 
describe it one must split the rapidities {uj} into those 
carried by the scalars, fermions and antifermions, which 
we denote as { 5 i},{ui} and jui}, respectively. The re¬ 
maining rapidities - those of the gluons - do not enter 
the matrix part since gluons do not carry R-charge. Then 

TT dwj r - 1-4 dwf 

- I -!- 2-k J 2tt 

i—1 i—1 

g(wi)g(w^)g(w3) 

/(wl, w2)/(w2, w3)/(wl, v)/(w2, s)/(w3, v) ’ 

where all integrals are over the real axis and where we 
use the short-hand notation 

ff(w) = J|g(Wi - Wj) , /(W,v) = Y[f{w^ - Vj) , 

i<j hi 

with f{x) = a;^ -b 1/4 and g{x) = x^{x‘^ + 1) to describe 
the fully factorized matrix part integrand, also depicted 
in figure 0 The numbers of auxiliary roots Rii, 2,3 are 
given by the solution to 

N^- 2 Ki+ K2 =Srb, 3 . 

+ Ki- 2 K 2 + K3= drb,2 , 

+ K2- 2K3 = Sr,,l . 
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All integrals in the matrix part representation can be 
straightforwardly performed by residues. For an NMHV 
hexagon with = 2, = 3 and = 0, for 

instance, we readily reproduce (§ in this way. 

This is it. All these ingredients can now be straight¬ 
forwardly plugged in Q to yield a fully non-perturbative 
representation of any six-point scattering amplitudes in 
planar Af = A SYM theory at any value of the coupling. 

As usual with bootstrap based approaches, our result 
lies somewhere between a physics conjecture and a math¬ 
ematical theorem. Any efforts to rigorously establish Q 
or to yield further evidence for it would be most valuable. 

In the meantime, it is a fascinating problem to inves¬ 
tigate the POPE sum in detail and look for physically 
interesting regimes where it simplifies. 

One such limit is weak coupling. At weak coupling al¬ 
most all integrals can be trivially done. More precisely, 
all fermion integrals over the so-called small fermion do¬ 
main can be performed straightforwardly leaving one 
with the integration over the rapidities of the others exci¬ 
tations - be them large fermions, scalars or gluons. The 
number of such integrations grows very slowly in per¬ 
turbation theory. For MHV say, at I loops, we have at 
most non-trivial integrations to perform! This is 
extremely economical as compared to any perturbative 
expansion and is one of the reasons why the POPE ap¬ 
proach is so convenient as a boundary data provider for 
the so-called hexagon bootstrap program [HHTH] . Re¬ 
producing the full kinematical dependence of the Tloop 
amplitude, away from the collinear regime, still requires 
summing over infinite families of flux tube excitations 
(arising mostly from integrating out any number of small 
fermions as discussed above). This calls for the develop¬ 
ment of efficient methods for both performing the OPE 
integrals and re-summing them, as tackled recently in 
[m [20] for the so called gluonic double scaling limit [7| . 

Another very interesting regime is strong coupling. 
Equipped with the POPE representation 0 one might 
want to derive the full strong coupling classical result mi- 
[25] . and address the question of computing sub-leading 
quantum corrections, following analyses of [6| [24H27] for 
amplitudes and similar problems. 

It is an equally fascinating problem to try to access 
other physically relevant kinematical regimes starting 
from the OPE sum. One such regime is the Regge limit 
[28{l3T] tackled from the POPE lens in [32l - l34| . It is also 
very interesting to establish a relation to the Q-bar ap¬ 
proach [23 [25] as initiated in m- 

Finally, we look forward to Monte Carlo simulating 
the flux tube gas Q on a computer and plot a six gluon 
amplitude all the way from weak to strong coupling [38j . 

Rare are the examples of CFT correlators for which 
the complete OPE series was found. Excitingly, thanks 
to integrablity, this turned out to be possible here. 
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